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Abstract 

We construct an explicit reversible symplectic integrator for the 
planar 3-body problem with zero angular momentum. We start with 
a Hamiltonian of the planar 3-body problem that is globally regu- 
larised and fully symmetry reduced. This Hamiltonian is a sum of 
10 polynomials each of which can be integrated exactly, and hence a 
symplectic integrator is constructed. The performance of the integra- 
tor is examined with three numerical examples: The figure eight, the 
Pythagorean orbit, and a periodic collision orbit. 

Keywords: 3-body problem; explicit symplectic integration; sym- 
metry reduction; Hamiltonian system. 



1 Introduction 

It is well known that the flow of a Hamiltonian H of the form H = 
T{p) + V{q) with conjugate variables q and p can be approximated by split- 
ting it into the integrable flow of T{p) and the integrable flow \l/y of V{q) 
and observing that = ^^o^*, + 0(t2), see, e.g. |2l El [El HO] . Thus a first 
order explicit symplectic integrator is obtained, and higher order methods 
can be constructed along similar lines [22] ■ In an integrable and separable 
Hamiltonian system of the form H = Hi{qi,pi) + H2{q2,P2) such splitting 
gives the exact identity = o "^^j^. If instead the Hamiltonian is a 
product H = Hi{qi,pi)H2{q2,P2) again the system is integrable with inte- 
grals Hi and H2 and the flow can be written as 
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In the superscript tHi denotes multiplication of t by the (constant) value of 
Hi. A monomial Hamiltonian is a special case that has the same structure. 
There are obvious generalisations to more degrees of freedom. Hence any 
polynomial Hamiltonian is a sum of integrable monomial Hamiltonians, and 
thus a splitting integrator can be constructed. Symplectic integration of 
polynomial Hamiltonians has been discussed in [161 SI El ID [IS] . 

In this paper we apply these methods to the polynomial Hamiltonian of 
the globally regularised and symmetry reduced 3-body problem at angular 
momentum zero. For a review of numerical and regularisation methods in 
the n-body problem we refer to p]. It is well known that binary collisions in 
the 3-body problem can be regularised. Regularisation consists of a canonical 
transformations which essentially extracts a square root near collision, and 
of a scaling of time so that the approach to the collision is slowed down. The 
classical simultaneous regularisation of the (spatial) 3-body problem is due to 
Heggie [7]. This increases the dimension of phase space from 18 to 24. Instead 
we would like to decrease the dimension of phase space by using reduction 
at the same time as regularisation. The simultaneous regularisation of the 
planar 3-body problem is due to Lemaitre [TT], and we use a version due 
to Waldvogel [21]. This is a symmetric simultaneous regularisation of the 
symmetry reduced planar 3-body problem and has the smallest possible 6- 
dimensional phase space. The resulting Hamiltonian is a polynomial of up to 
degree 6 in the canonical variables. A modern extension of these regularising 
transformations has recently been given by Montgomery and Moeckel [T3] . 
however, their Hamiltonians are not polynomial but rational. 

Our paper applies the methods for construction of an explicit symplectic 
integrator to Waldvogel's Hamiltonian with angular momentum zero. We 
also describe how a similar integrator could be constructed for Heggie's 
Hamiltonian, which works for non-zero angular momentum and in the spatial 
problem. 

2 The 3-body Hamiltonian 

The classical 3-body problem has long been studied, but still many open 
questions regarding its dynamics remain. For many questions, e.g. the study 
of relative periodic orbits, it is useful to reduce by translational and rotational 
symmetries, so that the absolute rotation of an orbit can be separated from 
shape dynamics in the centre of mass frame. Moreover, to study collision or 
near-collision orbits it is essential to perform (global) regularisation of the 
binary collisions. Following Waldvogel [21] we are going to do both. 

If the position and momentum of mass m^, for j = 1,2,3, are given by 
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Figure 1: Coordinates of a triangle in the plane with centre of mass at the 
origin O. 

complex Cartesian coordinates Xj and Pj respectively, we can transform into 
symmetry-reduced coordinates such that 

Xi - Xk = aje'f'^ , 

where aj = \Xi — Xk\ is the length of the triangle's side opposite to nij, (pj 
is the angle of that side in the original coordinate system (in the direction of 
nik to mi), as illustrated in figure [T| and j, k and / are cyclic permutations 
of (1,2,3). 

This reduction results in coordinates aj and = |(0i + 02 + 03), which 
represents the orientation angle of the triangle with respect to the original 
choice of Cartesian coordinates, and corresponding canonical momenta pj 
and p^. The Hamiltonian rewritten in these coordinates is independent of 
0, so p0 is a constant of motion. Hamilton's equations for {aj,(j),pj,p^) give 
the reduced dynamics, including a differential equation for which may be 
integrated along to be able to recover the unreduced position of the triangle. 

The globally regularising transformation, illustrated in figure |2} goes from 
symmetry-reduced to regularised coordinates, simultaneously regularising all 
the binary collisions. Define aj for j = 1, 2, 3 such that aj = a\ + af. In this 
way a| is the distance from to the point where the incircle of the triangle 
touches the sides adjacent to mj. 
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Figure 2: Physical significance of tlie regularised coordinates aj. 

The space of coordinates (oj) is the space of all triangles, not accounting 
for orientation. Orientation is taken to be positive if, going clockwise around 
the triangle, the masses are encountered in a cyclic permutation of (1,2,3) 
or negative otherwise. The space of all possible oriented triangles is called 
the shape space, and the space of (aj) is a four-fold covering of this space, 
in which the sign of the product aiQ;2tt3 determines the orientation of the 
triangle. Thus the triangle formed by {ai,a2,a3) is the same as the ones 
formed by (ai, —ck2, —eta), (— ai, a2, —0^3) and {—cki, —a2, eta). 

Canonically conjugate momenta ttj are introduced using a generating 
function. Finally the time scaling 



together with Poincare's trick to make this Hamiltonian yields the regularised 
and symmetry reduced polynomial Hamiltonian 



where H is the original Hamiltonian written in the new coordinates and 
h = H{oio, TTo) is the energy corresponding to the initial conditions (ao, ttq), 
so only those solutions for which K = are physically meaningful. 

The Hamiltonian of the zero-angular momentum 3-body problem in reg- 
ularised coordinates is 



dt 

— = 010,203 



(1) 



K ^ {H - h)aia2a3, 



K ^ Kq — haia2a3, 



(2) 



where 




(3) 
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in which 



TT = [TTi TT2 TTs j 

B{cx) = ^3 A2 5i , 
\B2 Bi As J 

where 

Aj = ^a' + ^af + ^al 
ruj ruk mi 

Bj — ——a^ai and 

2 9 2 2 



The new Hamiltonian is a polynomial in ol and tt, and thus Hamilton's 
equations of motion for this system can be integrated with an explicit sym- 
plectic integrator obtained by splitting into monomials. As we are going to 
show in the next section it is more efficient to split into certain polynomials 
whose flow can be exactly solved. 



3 Construction of the Symplectic Integrator 

The basic building block of the integrator are the exact solutions for mono- 
mial Hamiltonians i^^n = q"^p^ in one degree of freedom. The flow of this 
Hamiltonian for m 7^ n is 

= where ^ = (l + (n - m)q^-^p^-H) — (4) 

while for m = n it is 

^l{q,p) = {qP,p/Pl where ^ = exp(m(gp)"^-ii) . (5) 

For the Hamiltonian we are studying the cases that occur are n — m — 1, 
n — m — 2, and m — 3, n — 1. We also recall that if the Hamiltonian 
is a function of positions or momenta only (with any number of degrees of 
freedom) the flows are 

^r(p) (9, P) = {q + ^pT, p) , and {q,p) = {q,P - ^qU) . (6) 
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3.1 Integrable Polynomial Hamiltonians 

The basic building blocks just mentioned are now combined to form integra- 
tors for the terms that appear in the Hamiltonian K. The main observation 
is that if the Hamiltonian is a product of factors that depend on disjoint 
groups of degrees of freedom, then each factor is a constant of motion. Each 
of the factors in our case is either depending on momenta or positions only 
(denoted by T{p) or U{q)) or it is a single monomial in one degree of free- 
dom (denoted by Hmn) or a sum of monomials of disjoint degrees of freedom 
(denoted by G). 

We now list the cases that are relevant in our case (recall that each of the 
factors depends on disjoint groups of degrees of freedom): 

Ha = THnm, = V't^"™ O ^^'^ (7a) 

H, = TV, ^, = V'^^o7/>f (7b) 
H, = GHnm, ^c = ^f^o^'^- (7c) 

where G is a Hamiltonian which is the sum of Hamiltonians depending on 
disjoint degrees of freedom G = Hi{qi,pi) + -^2(92,^2) and thus ipG = i'Hi ° 
ipH2- Note that all these formulas are exact, and that the order of composition 
is irrelevant since the flows commute and the individual factors are constants 
of motion. 

3.2 Splitting 

Let us now explain how to split K ^ into such terms. It is a polynomial 
Hamiltonian of degree 6 in ck and tt with 34 monomials. There are 13 
monomials, dependent only on a., of degrees 6 and 4, which may be treated 
as a single stage. The remaining 21 terms may be grouped such that only 
9 more stages are necessary to approximate the flow of the full Hamiltonian 
to first order in the time step in 10 stages. Let K = ^2^=0 "where we set 
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Mj = rukTrii and = + Then the sphtting is 



- Mja^ - (J2 Mj) (E - haia2a3 = Cq 

I {N2al + Nsal) ol\ti\ = | C1.23 

i {Nsal + N,al) aj^l = | Cs.ai ^2,2 

(iVia? + N^al) ajTil = | 6-3,12 63,3 

^3«3 + 4«3«? + V2 = (8) 

'N,a\ + ^alal + N2at)Til = | 

;^a2vr2 + i^oisTT^] al-Ki = -| 67,23 67,1 

:;^(^3'^3 + ^(^I'^i) "2^2 = -| 63,31 63,2 

^"1^1 + ^;!T'^27r2j afTTs = -| 69,12 69,3, 

where each subindexed function 6j is a constant of motion in its associated 
Hamiltonian. 

There are clearly four groups in equation (|8]), which we shall enumerate 
0: {0}, 1: {1,2,3}, 2: {4,5,6} and 3: {7,8,9}. Hq depends on coordinates 
only, so can be integrated by Group 1 can be integrated by (7b), group 




2 can be integrated by (7a), and finally group 3 can be integrated by (7c) 



where 6 is a sum of H,mm Hamiltonians. 



3.3 Higher order methods 

An important ingredient in constructing higher order reversible methods is 
the adjoint (0*)* of a method 0* which is defined to be If 0* = 

■0^ o o ■ ■ ■ o V'n each -0* is self-adjoint, then the adjoint is obtained by 
reversing the order of composition. This follows from the definition of the 
adjoint: 

(0*)* = (0-*)-^ 

In our case the self-adjointness of the individual steps ipj follows from the 
fact that they are exact solution of Hamilton's equations. 

Channell & Neri [2] offer a basic derivation of a reversible, symplectic 
map that is accurate to second order in the time step. When the splitting is 
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of the form H = T{p) + U (q) this leads to the symplectic leapfrog integrator, 
by composing symplectic Euler with its adjoint. 

This construction also applies to the more complicated case with a first 
order integrator composed of 10 self-adjoint maps as in our case. Given 
(j)'^ = ipl o o ■ ■ ■ o ip^ as above a reversible second order method is found as 

02 = 4 ° = O . . . O ^|_, O tPl^ O O . . . O . (9) 

Yoshida's 1990 paper [22J gives a general method by which one may obtain 
integrators of arbitrary even order, if only one has, to start with, a reversible 
even-order integrator 02 such as the symplectic leapfrog — or, more generally, 
symplectic midpoint. One can compose (pl to obtain a fourth order integrator 
04, and compose these to obtain 0g and so on. In general, given 02n5 

02n+2 = 02n 02n 02n ) (10) 
2l/(2n+l) I 

where we define zq = - 2i/(2(2«+i)) ' = 2 - 2V(2»+i) *° ^^i^^^ ^tep 
size of the lower order method. 



3.4 Finite time blowup 

It must be noted that the solution of the Hamiltonian H = q^p^ given in Q 
can (for n ^ m) reach infinity in finite time. This occurs when 

1 + (n - m) g^-i p^^^^H = 0. (11) 

This obviously makes step sizes comparable to this threshold risky when 
this form of solution is used in the integrator. This singularity could be 
reached if the denominator is negative and large during forward timesteps, or 
if the denominator is positive and large for "backward" timesteps (as during 
the middle stage of Yoshida's trick). This possibility arises in equation ^ in 
the group 4 of the splitting ([s]), which have terms of the form a^vrj. It may 
appear that this finite time blowup is an artefact of the integrator. However, 
after the time scaling the Hamiltonian K does have finite time blow up when 
particles escape to infinity. In this light it seems less unexpected that a stage 
of the corresponding symplectic integrator shows the same behaviour. 

Blanes [1] provides a means by which to avoid such singularities, by way 
of rewriting the polynomial in terms of sums of binomials in the coordinates 
and momenta and finding coefficients such that the two expressions are equal. 
We did apply this to the Hamiltonian H31 in our problem, and found a way 
to replace this with a Cremona map. However, it turned out that the overall 
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error of the method was worse than without this modification. Our method is 
more expensive, since it needs to compute rational powers, but this additional 
cost is worth it. 



3.5 Other polynomial globally regularised Hamiltoni- 
ans 

Our main concern in this paper is the zero-angular momentum reduced and 
regularised planar 3-body problem, which has 3 degrees of freedom. Other 
well known globally regularised polynomial Hamiltonians are due to |2D] for 
the planar 3-body problem and to [7J for the spatial 3-body problem. These 
Hamiltonians are not fully symmetry reduced and have 4 degrees of freedom 
(planar arbitrary angular momentum) and 12 degrees of freedom (spatial 
arbitrary angular momentum). 

Reggie's simultaneously regularised Hamiltonian for the spatial 3-body 
problem [7j in canonical variables Qji and conjugate Pjj, j = 1,2,3, i = 
1, . . . , 4 has the form 

H = Hq + H4 + ifs + Hq + 
Hq = —haia2a3 — mjUikajak, 

H3+^ = l^H', 1 = 1,2,3 
4 ^-^ rrij 

where fijk = mjUik/imj + ruk), and aj = Y.t=iQ]i^ biP = ZlLi -P/r 1^ 
addition Pi = {Pn, P12, P13, Pu)^ and Ai is the KS-matrix [9J of the form 



The terms Hq and -^4,5,6 are analogous to the previous ones. The terms in 
H _i can be split into 9 terms of the form ajf^gi where the functions fk and 
Qi only depend on the degrees of freedom k and I, respectively. These terms 
are somewhat similar to the Hamiltonians -^1,2,3 in Waldvogel's case. Thus 
the Hamiltonian can be split into 13 polynomials of degree up to 6 each of 
which is integrable. 

When setting Qji and Pji with i = 3, 4 equal to zero Heggie's Hamiltonian 
describes a planar problem. However, this still has 6 degrees of freedom. We 
can reduce the number of degrees of freedom to 4 by instead using Waldvogel's 
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Average single timestep error vs timestep size 




Figure 3: Performance of the integrator at 2nd, 4th, 6th and 8th orders aver- 
aged over a large number of points in phase space at varying time step sizes. 
Evaluation time for a fixed number of steps is indicated. These computations 
were performed in Matlab R2012a (win64) on a quad core Intel Core 17 Q 
820 processor at 1.73GHz. 

Hamiltonian [201 E]- This Hamiltonian is a polynomial of degree 12 and can 
be split into 15 terms in a way similar to the two cases discussed above. 

4 Numerical examples 

In this section we will show some numerical results achieved using our integra- 
tor in a selection of orbits ranging from far from collision to close encounters 
to a collision orbit. 

Figure |3] shows the single-timestep energy error (averaged over a large 
volume of phase space comprising 40649 initial conditions) of the integrator 
at different orders of accuracy and different time steps. It also shows the 
amount of time and total number of steps taken in reducing step sizes until 
the error falls below 10~^^. From this we can choose an order that is most 
suitable to our needs, optimising speed and accuracy. 

From the previous section, the number of evaluations of the base second 
order method at order 2n is 3""^. In order to integrate a given interval of 
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time with average single timestep error below 10^^^, the data in figure [s] 
suggest choices of 10"^, 10~^, 10~^ and 10~^-^ ^ 3 x 10"^ for orders 2, 4, 6 
and 8 respectively. To integrate for a time interval of 1 with those step sizes, 
at second order we must evaluate the base method 10^ times, at fourth order 
3 X 10^ times, at sixth order 3^ x 10^ = 9 x 10^ times, and at eighth order 
3^ X I X 10^ = 9 X 10^ times. Clearly the gains in accuracy start to become 
outweighed by the burden of complexity by the eighth order, so we choose 
the sixth order algorithm to integrate our examples. 

Consider the figure-8 choreography, discovered by Christopher Moore 
|14j . proved to exist by Alain Chenciner and Richard Montgomery [3] and 
explored by Carles Simo [UlITTj. We choose initial conditions 

ao = (0,1.134522804969261,1.134522804969261)^ 

TTo = (1.506773685132772,0.694233777317562,-0.694233777317562)^ 

in regularised coordinates, with h = —1 and equal unit masses. Integration 
is using a 6th order algorithm with time steps Ar = 0.001. In scaled time, 
the figure — 8 has a period of 2.221813718; in physical time its period is 
9.2371333. The trajectory in regularised coordinates is shown in figures |4a] 
and |b] and the energy error over 25 orbits is shown in figure |4cj Figure [5] 
shows the trajectory of this orbit in the 3-dimensional space {aj). Note that 
each crossing of a plane aj = corresponds to a syzygy with rrij in the middle 
of the configuration. 

Next we look at the Pythagorean orbit [I9j for mi = 3, m2 = 4, = 5, 
with initial conditions as given in [S], which, in regularised coordinates, are 

ao = (1, VS, V2f 
TTo = (0,0,0)^. 

This orbit has a close encounter between masses 1 and 3 at around t = 15.8 in 
physical time (about r = 1.52 in scaled time). In order to preserve accuracy 
in the "energy" for this trajectory, we integrate at order 6 with time steps 
Ar = 0.0001, owing to the larger masses and greater distances between 
bodies on average. Waldvogel's analysis regularises the system, albeit slightly 
differently, and his integration is not symplectic. The final motions of this 
orbit compare well with other studies; plotting the orbit in physical space 
produces results indistinguishable from [19], [5]. 

The regularisation of the 3 — body problem allows our integrator to cope 
well when the distances between any two masses are small. The result of 
the time scaling is that the regularised system has a finite time blowup for 
any escape orbit. If one continues to integrate the Pythagorean orbit past 
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(a) Figure-8 choreography in regu- (b) Regularised momenta of figure-8 
larised coordinates integrated at 6th choreography integrated at 6th order, 
order, At = 0.001. At = 0.001. 



Error between forward and backward integrations 





(d) Two-way integration error at 6th 
(c) Energy error integrated at 6th or- order, showing the error when inte- 
der for 25 periods with 200 time steps gration is reversed back to r = dur- 
per period. ing the 25-period integration. 



Figure 4: The figure-8 choreography. 



T = 8.105 with At = 0.0001, the error in the energy grows exponentially and 
the results become inaccurate. 

Finally, we show results in a periodic collision orbit, discovered during a 
search for periodic orbits in the reduced space, with initial conditions 

cxo = (0,0.717162073833634,1.683647749751810)^ 

TTo = (1.762174970761679,0.177158588505747,-0.401743282150556)^, 

equal unit masses and h = —1, integrated at 6th order with At = 0.001. This 
orbit has collisions between masses 1 and 2, as can be seen by ai = 0^2 = 



at r = 1.9362 and r = 5.062 in figure 7a This orbit is shown in figures ^ 



and^ Its scaled period is 6.2520511, corresponding to a physical period of 
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Trajectory in regularised coordinates 




Figure 5: Trajectory of figure — 8 choreography in a-space, lying nearly on a 
plane. Colour gradient represents the moment of intertia (lighter is higher). 



29.6117209. Note in figure |9] that the collisions happen on the a^-axis, when 
0(1 = 0(2 = 0. 

Because Hamiltonian systems are time-reversible, it is desirable to have an 
integrator with the same property. The second order map is constructed 
as such, and Yoshida's formula for higher order integrators constructs them 
to be reversible as well. That means that 02n 02n — roundoff error. 

Figures 4d, 6d and 7d show how closely this integrator returns to its initial 
condition after a certain number of time steps in one direction, followed by 
the same number of iterations with a negative time step. 

By this measure, the integrator has the most trouble with Pythagorean 
orbit, which clearly shows signs that it exists within a chaotic region of phase 
space by the exponential growth of error. However; energy is well preserved 
in this and the other cases. 
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2 3 4 5 



(a) The Pythagorean orbit in regu- (b) Regularised momenta of the 
larised coordinates. Pythagorean orbit. 




1 2 3 4 5 6 7 



(c) Energy error. 



Error bebween forward and backward integrations 




(d) Two-way integration error. 



Figure 6: The Pythagorean orbit integrated at 6-th order with time steps of 
0.0001 up to r = 8.105 in scaled time. 
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(a) Periodic collision orbit in regu- (b) Regularised momenta of periodic 
larised coordinates. collision orbit. 



Energy Error between forward and backward integrations 




(c) Energy error for 5 periods of the (d) Two-way integration error over 5 
periodic collision orbit. periods of the periodic collision orbit. 

Figure 7: A periodic collision orbit integrated at 6-th order with time steps 
of 0.001 in scaled time. 
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Cartesian reconstruction 
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Figure 8: Reconstruction of Cartesian trajectories from regularised integra- 
tion for a periodic collision orbit. 



Trajectory in regularised coordinates 




Figure 9: Trajectory of periodic collision orbit in a-space. Colour gradient 
represents the moment of intertia of the configuration (lighter is higher). 
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5 Conclusion 



We have constructed a symplectic integrator for the reduced and regularised 
planar 3-body problem at zero angular momentum. The method works well, 
but it is not very efficient, because each (ffist order) time step involves the 
computation of 10 individual maps. Our interests is the computation of 
relative periodic orbits including collision orbits, and for this task the method 
is appropriate. The detailed results about relative periodic orbits and their 
geometric phase will be reported in a forthcoming paper. 



A Integrator stages 



Subsection described how to integrate a monomial Hamiltonian, and sub- 
section 3.2 described the splitting of equation ^ into a minimal number of 
solvable parts and those solutions. Here we use those solutions to build an 
explicit ffist order symplectic composition method for ([2]). 

Let the timestep be At, let /ij = {rrik + m^), let the values of the system 
before and after one timestep respectively be Zq = {ai^o, . . . , tts o) 

^'^ and intermediate steps be = . . . , 7r3^.j_i)-^. 



. . . ,7r3,i, 



''^ and Zi 



Now 



^1 



"1,0 
"2,0 
"3,0 

TTi^o + 2ai,o ((2"?,o + «i,o) (^^1,0 + Mi) + miai,o/iiai,o) Ar 
7r2,o + 2^2,0 ((2"|o + a2,o) {ha2,o + M2) + m2"2,oAt2a2,o) Ar 
\tt3^o + 2a3fl ((2a|o + 03,0) {has^o + M3) + m3Q;3,o/i3a3,o) At J 

/ exp (1 (A^2"2,.i + ^3"I,.i) "i,.i TTi,.! Ar) ^ 

"2,.l 
"3,.l 

TTi,.! exp (-| (A^2"i,.i + A^3"i,.i) "i,.i Ar) 
7r2,.i -\N2 al l '^2,.iAr 
V 7r3,.i - I A/'s al l "3,.iAr / 

/ «1,.2 \ 

"2,.2 exp {\ (A^3"3,.2 + ^l"l,.2) "2,.2 7r2,.2 Ar) 
"3,.2 

7ri,.2 -\Ni a| 2 ^2,.2 "i,.2Ar 
vr2,.2 exp {-\ (A^3"3,.2 + Nial J) a2,.2 7r2,.2 Ar) 



7r3,.2 



i N3 "1.2 



TT. 



2, .2 "3, .2 



Ar 
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/ «1,.3 \ 

"2,.3 

q;3,.3 exp (l (iViaf 3 + N2al^) as^.s 713^,3 At) 
7ri,.3 - i iVi 7r| .3 ai,.3Ar 
7r2,.3 - J ai,.3 7r|,.3 a2,.3AT 
\7r3,.3 exp {-I {Nial -i + Ar2Q;| 3) q;3,.3 7r3,.3 At) J 

^ai,.4 + 7ri,.4 (1 (iV2a|.4 + iV3tt3,.4) + i«i,.4«i4) 

q;2,.4 
a3,.4 
7ri,.4 

7r2,.4 + Itt? .4 fA^2a|.4 + ;^Q;2,.4Q;i,.4) At 
\^ 7r3,.4 + ^TT^ 4 [N^al^ + ;^«3,.4a;i,.4j At J 

( «1,.5 \ 

«2,.5 + 7r2,.5 (i (iV34.5 + ^i^.s) + s^^^als^ls) At 

a3,.5 

7ri,.5 + ^TtI .5 (^1«?,.5 + ;4"l.-5"i,-5) A^ 
7''2,.5 

^ 7r3,.5 + |7r| .5 (^^3ai,.5 + ^^"s.-s^ls) A^ y 

/ «1,.6 \ 

«2,.6 

a3,.6 + 7r3,.6 (i (iVi<.6 + ^2a2,.6) + ai^ale^le) 
7ri,.6 + ivrl .6 f^^i"?,.6 + ^"I'-e^Lfi) 

V 7r3,.6 / 

'^ai,.7 (1 + \ (4tt2,.77r2,.7 + 4a3,.77r3,.7) a?,.7AT) 
a2,.7exp f-4^af 77ri,.7AT 
a3,.7exp(^-4^a;3.77ri,.7AT 

^lv7 (1 + ^ (4«2,.77r2,.7 + 4«3,.77r3,.7) a?,.7A^) 

7r2,.7exp f^^a? 7^i,.7At 
7r3,.7exp (^i^a^ 7^^^^^^^ 
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/ ai,.8exp (^-^a| 87r2,.8Ar) \ 

Qi2,.8 (l + I (^Q;3,.87r3,.8 + ^^ai^sTTi^g) al.gAr) 
Q;3,.8exp(^-^Q;|87r2,.8AT) 
^i.sexp r^^^ai 87r2,.8AT) 

3 

^2,.8 (l + i (;i^7a3,.87r3,.8 + ^^ai„s7ii^,s) ai,.8AT) ' 
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^3,.8exp (^j^Q;|.87r2,.8AT) 

Q;i,.9exp f-^af 9^3,.9At 
a2,.9exp (^-^a| 97r3,.9Ar 

«3,.9 (l + I (4«i.-97ri,.9 + ^Q;2,.97r2,.9) "l^gAr) 

7r2,.9exp (^^af 97r3,.9AT 
\^ 7r3,.9 (l + I (4"i,.97ri,.9 + ;^a2,.97r2,.9) al.gAr) ' J 
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